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SUMMARY 


The  purpose  of  the  reported  research  is  to  study  the 
dynamics  and  control  of  a  class  of  large  antenna/reflector 
systems  in  orbit  which  are  also  partially  stabilized  using  a 
tether-connected  subsatellite.  The  initial  focus  has  been 
in  the  development  of  the  system's  equations  of  motion 
linearized  about  the  equilibrium  position  where  the 
reflector's  (shell's)  symmetry  axis  nominally  follows  the 
local  vertical.  The  shell  roll,  yaw,  tether  out-of-plane 
swing  motion  and  out-of-plane  elastic  vibrations  are 
decoupled  from  the  shell  and  tether  in-plane  pitch  motions 
and  in-plane  elastic  vibrations.  It  is  proved  that  the  in¬ 
plane  motion  of  the  system  could  be  asymptotically  stable 
based  on  Rupp's  tether  tension  control  law  based  only  on 
length  and  length  rate  information.  However,  the  transient 
responses  can  be  improved  significantly  (especially  for 
damping  of  the  tether  and  shell  pitch  motion)  by  using  an 
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by  different  investigators  has  been  completed.  Recommen¬ 
dations  are  made  concerning  the  suitability  of  the  various 
control  laws  for  use  with  the  orbiting  tethered  reflector 
system . 

Finally  a  preliminary  model  of  the  nonlinear  dynamics  of  the 
tethered  antenna/reflector  system  in  orbit  has  been  obtained 
based  on  Lagrangian  techniques.  It  is  seen  that,  unlike  the 
situation  for  the  system  linearized  about  the  nominal 
stationkeeping  motion,  the  in-plane  and  out-of-plane  motions 
are  coupled  through  second  order,  and  nonlinear  coupling 
terms  also  depend  on  tether  line  swing  motions  and  tether 
vibrations.  For  this  preliminary  model  the  shell  is 
considered  to  be  a  rigid  structure. 
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1 .  INTRODUCTION 


1.1  Feasibility  of  concept  based  on  Existing  Work 

Since  the  early  1970 ' s  a  number  of  very  large  space  antennas 
have  been  proposed  for  power  transmission,  astronomical 
research  and  communications.  The  gravity  stabilized 
configuration  is  particularly  suited  for  a  very  large 
flexible  structure  to  alleviate  the  problems  associated  with 
the  active  attitude  control  of  very  large  structures.  The 
structural  feasibility  of  a  very  large  Earth  oriented 
antenna,  where  the  flexible  reflector  contour  is  maintained 
by  adjusting  the  length  of  connecting  tethers  between  the 
reflector  and  feed  panels,  has  been  discussed. t 1 1  In  this 
paper  the  stress  analysis  of  the  tethered  antenna  was  given. 
The  analysis  of  the  dynamics  and  control  of  the  orbiting 
flexible  shallow  spherical  shell  and  various  tether 
connected  systems  in  space  have  been  performed.  Bainum  and 
Kumart^J  have  investigated  the  dynamics  of  an  orbiting 
flexible  shallow  spherical  shell  with  a  dumbell  connected  to 
the  shell  at  its  apex  by  a  spring-loaded  double-gimball 
joint  to  provide  the  favorable  composite  moment  of  inertia 
distribution.  Also,  Bainum  and  Reddy^J  have  investigated 
the  shape  and  orientation  control  of  this  shell  antenna  by 
including  some  additional  active  control  elements. 

Numerical  results  verify  that  a  significant  savings  in  fuel 
consumption  can  be  realized  by  using  the  hybrid  shell- 
dumbell  system  together  with  the  (active)  point  actuators. 
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The  purpose  of  the  proposed  research  is  to  study  the 
dynamics  and  control  of  a  class  of  large  antenna/reflector 
orbiting  structures  which  include  an  articulated  tether 
connected  supporting  structure  to  provide  the  favorabje 
moment  of  inertia  distribution  for  over-all  gravitational 
stabilization  together  with  some  active  actuators.  There 
are  two  possible  proposed  subsystems  which  could  provide  the 
connection  between  the  tether  and  the  shell  reflector;  one 
involves  a  spring-loaded  doubled-gimballed  joint  connected 
to  the  shell's  apex  and  through  which  the  tether  is 
deployed/retrieved  (Fig.  1);  the  second  contains  a  joint  at 
the  end  of  a  rigid  boom  which  is  attached  to  the  shell's 
apex  (Fig.  2).  Through  the  end  joint  the  tether  would  be 
deployed  or  retrieved.  The  tether  tension  could  be  used  for 
producing  restoring  torques  on  the  shell,  with  natural 
damping  provided  in  the  joint  assembly.  For  the  first  phase 
of  the  study  reported  here  the  second  subsystem  has  been 
taken  as  the  basis  for  the  system  model  due  to  the 
relatively  simpler  implementation  as  compared  with  the 
double-gimballed  joint  in  the  first  subsystem. 

1.2  Relevance  to  SDI 

Associated  with  the  capability  to  orient  a  large  flexible 
antenna/reflector  type  of  device  accurately  while  at  the 
same  time  maintain  the  surface  shape  tc  within  centimeters 
or  even  millimeters  are  many  applications  in  both  the 
military  and  rivilian  fields.  For  example,  high  energy 
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Fig.  1.  Subsystem  A  -  Tether  Deployed  from  Apex  of  Reflector 
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Subsystem  B  -  Tether  Deployed  from  the  End  of  a 
Rigid  Boom  Connected  to  the  Reflector 
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beams  can  be  generated  by  a  power  source  and  reflected  from 
specific  known  points  on  the  reflector  surface  to 
preselected  targets.  In  the  very  important  communications 
field,  such  an  antenna  surface  can  receive  multibeam 
communication  waves  from  electronic  feed  devices  and 
transmit  these  to  a  variety  of  small  mobile  receivers  to 
comprise  strategic  communication  links  during  early, 
critical  phases  of  an  attack  when  larger,  fixed  land-based 
antennas  would  be  far  more  vulnerable  to  observation/damage. 
Such  devices  could  also  be  employed  to  transmit  coded 
electronic  mail  rapidly  over  different  communication 
channels . 

1.3  Outline  of  the  Research  Reported 

The  second  chapter  focuses  on  the  development  of  the  linear 
system  equations  of  motion  for  an  orbiting  tethered  shallow 
spherical  shell  system  where  the  shell's  axis  of  symmetry 
nominally  follows  the  local  vertical.  The  Newton-Euler 
method  for  a  continuous  system  is  adopted  here.  The  second 
objective  is  to  develop  the  in-plane  and  out-of-plane 
stability  conditions  and  introduce  some  tension  control  laws 
for  in-plane  motion  control.  The  transient  responses  will 
be  compared  for  three  different  tension  control  laws  during 
typical  station  keeping  operations.  A  paper  based  on  these 
tasks  was  presented  at  the  Third  international  Conference  on 
Tethers  in  Space,  San  Francisco,  Masy  17-19,  1939  and  has 
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been  accepted  in  slightly  revised  format  for  publication  in 
The . Journal . .of . the ..  Astronaut  real. .  S  c.i  erLc.es.... 

The  following  chapter  describes  a  comprehensive  review  of 
the  steps  in  the  development  of  control  laws  for  Shuttle  or 
platform  connected  tethered  subsatellite  systems. 

Deployment,  stationkeeping,  and  retrieval  control  strategies 
are  reviewed  and  compared.  Finally,  recommendations  are  made 
suggesting  the  relative  suitability  of  the  different  control 
laws  for  adaptation  with  the  proposed  orbiting  tethered 
reflector  systems. 

Chapter  Four  concentrates  on  the  development  of  the 
nonlinear  equations  of  motion  for  the  tethered  reflector 
system  in  orbit  m  a  form  suitable  for  simulasting 
deployment  and/or  retrieval  maneuvers  based  on  some  of  the 
control  laws  described  in  Chapter  Three. 

Finally,  Chapter  Five  summarizes  some  concluding  statements 
and  follow-on  plans  for  the  continuation  of  this  general 
area  of  research. 
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2  DYNAMICS  AND  CONTROL  OF  A  TETHERED  ANTENNA/REFLECTOR 
IN  ORBIT 
2.1  Introduction 

Since  the  early  1970's  a  number  of  very  large  space  antennas  have  been 
proposed  for  power  transmission,  astronomical  research  and  communications.111  The 
gravity  stabilized  configuration  is  particularly  suited  for  very  large  flexible  systems 
to  alleviate  the  problems  associated  with  the  active  altitude  control  of  very  large 
structures.  Bainum  and  Kumar121  have  investigated  the  dynamics  of  an  orbiting 
flexible  shallow  spherical  shell  with  a  dumbbell  connected  to  the  shell  at  its  apex 
to  provide  the  favorable  composite  moment  of  inertia  distribution.  Also,  Bainum  and 
Reddy111  have  investigated  the  shape  and  orientation  control  of  this  shell  antenna 
by  including  some  additional  active  control  elements.  Meanwhile,  scores  of 
applications  of  tethers  in  space  have  been  proposed  and  analyzed  including  some 
space  platform-based  applications  of  the  tether  subsalcilitc  system. “MS1 

The  objective  of  the  present  paper  is,  first,  to  develop  a  system  mathematical 
model  of  a  class  of  large  anten.ia/rcflcctor  orbiting  structures  which  include  an 
articulated  tether-connected  supporting  structure  to  provide  the  favorable  moment 
of  inertia  distribution  for  over-all  gravitational  stabilization,  together  with  some 
active  actuators.  The  tether  would  be  connected  at  the  end  of  a  rigid  boom  which 
is  attached  to  the  shell's  apex  and  through  the  end  of  the  boom  the  tether  could 
be  deployed  or  retrieved  (Fig. 3).  The  tether  tension  could  be  used  for  producing 
restoring  torques  on  the  shell.  The  second  objective  is  to  develop  the  in-plane  and 
out-of-plane  stability  conditions  and  introduce  some  tension  control  laws  for  in-plane 
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motion  control.  The  transient  responses  for  the  three  different  tension  control  laws 
will  be  compared  during  typical  station  keeping  operations. 


2.2  Equations  of  Motion 

For  system  modelling  the  following  assumptions  were  made: 

1)  The  thickness  of  the  shell  is  small  as  compared  to  the  height  of  the  shell, 
and  the  ratio  of  the  height  to  the  base  radius  is  much  less  than  unity  (  condition 
for  shallowness  ). 

2)  The  elastic  deformations  perpendicular  to  the  symmetry  axis(  i.e.^x  axis  )  of 
the  shell  are  negligible  compared  with  the  deformations  parallel  to  the  symmetry 
axis,  i.e.,  only  transverse  vibrations  are  considered. 

3)  The  symmetry  axis  of  the  shell  is  nominally  along  the  local  vertical. 

4)  The  center  of  mass  of  the  system  is  moving  in  a  circular  orbit. 

5)  The  flexibility  of  the  boom  is  neglected. 

6)  The  subsatellite  is  to  be  considered  as  a  point  mass. 

The  shift  of  the  center  of  mass  of  the  system  will  be  considered.  In  order  to 
develop  a  general  model  for  the  tethered  shell  system  it  is  assumed  that  the 
massive,  flexible  tether  is  deploying  or  retrieving  a  subsatellitc  at  a  distance.  ,  from 
a  point  on  the  shell  which  is  offset  by  distance  h„.  h.,.  h,.  along  the  yaw,  pitch, 
and  roll  axes,  respectively,  from  the  center  of  mass  of  the  shell.  O,. 
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Santini1*',  Bainum  and  Kumar"'  have  developed  a  mathematical  formulation  for 
a  general  orbiting  flexible  body  based  on  the  Newton- Euler  method  and  continuum 
approach.  In  the  present  paper  this  method  will  be  extended  to  the  system 
composed  of  two  flexible  structures  (  the  shell  and  the  tethered  subsatellite  ). 

The  coordinate  systems  used  in  the  development  of  the  system  equations  of 
motion  are  shown  in  Fig.  4.  O.X«YaZa  is  an  orbit-fixed  reference  frame  centered  at 
the  center  of  the  mass  of  the  sheil,  0„  with  0,Xa  along  the  local  vertical  and  O„Y0 
along  the  orbit  normal  opposite  to  the  angular  velocity  vector.  O.X.Y.Z,  is  an 

undeformed  shell  reference  frame.  R„  where  O.X.,  O.Y,.  O.Z,  are  the  principal  axes 

of  the  sheil.  OXYZ  is  the  subsateliite-undeformed  tether  reference  frame.  R„  with 
OX  along  the  undeformed  tether  line,  where  O  is  the  point  from  which  the  tether 
is  deploying  or  retrieving.  The  coordinates  of  O  in  the  shell  frame,  R„  are  h„  h¥, 

K 

The  angles  v  ,  3  ,  b  are  the  yaw.  pitch  and  roil  angles  of  the  shell, 

respectively.  An  Euler  angle  rotation  sequence  of:  (I)  4>  .  (2)  9  .  and  (3)  is 

assumed  from  the  0,.\aY,Za  system  to  the  O.X.Y ,Z.  system. 

The  transformations  from  O,X0Y3Z0  to  O.X.Y.Z,  and  from  O.X.Y .Z.  to  OXYZ 
are  assumed  to  be  given  by 


* 

** 

'x. 

cbc9 

s$c'.{H-c$sds^ 

sosijj-cisdei) 

X, 

Y’ 

-scc9 

cbc'i-sbsQ  sty 

ces^+scisdc'i/ 

Y, 

1 

z. 

33 

-C3S0 

*».  W  W  ♦  ■ 

Lz3. 
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t  x  y  z  r 


T(a  y)  [  X.  Y,  Z.  ]r 


where 


-cyca 

sy 

sacy 

T(a,y)  - 

3yca 

cy 

-sosy 

-90 

0 

-co 

where 

c  ♦  cosine  (  ), 

s  * 

sine  (  ) 

(2) 


Consider  an  elemental  mass,  dm,  whose  instaneous  position  vector  from  the 
center  of  the  shell,  0„  is  r  (Fig.4).  The  equation  of  motion  for  dm  can  be  written 
as17’'1" 


a  dm  31  L(q)  +  f  dm  +  c  dm  (3) 

where  a  =»  inertial  acceleration  of  dm 

q  -  elastic  displacement  vector  of  dm 
L(q) «  elastic  forces  acting  on  dm 
f  =■  gravitational  force  per  unit  mass 
e  *■  external  forces  acting  per  unit  mass 
The  gravity  force  in  the  shell  frame,  R„  is  given  by:’ltn 

f  -  f0  +  M  r  (4) 

where  f,  is  the  gravity  force  at  O,  expressed  in  the  I'm  me.  R„.  and 


M 


3c2-ic2e-l  -3sa>cPc2e  3cpc9s8 
33<Dcbe2e  332‘5c26-1  -3sdcSsd 
3cbs9c9  -3soc9s0  3s29-1 


(5) 


where  •-»  is  the  orbital  angular  velocity. 
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The  vector  equation,  (3),  can  be  written  in  the  frame,  R„  as 
'*  •  • 

CaQ-f 0+r  +2o«r+5ac(aacx)+SxT-te]  dm-L(q)  -^dm  ■  0  (63 

•  *  * 

where  r,  r  are  the  velocity  and  acceleration  of  dm.  respectively,  as  seen  from 
the  frame,  R„  and  <J  is  the  angular  velocity  of  the  frame  R,. 


0) 

X 

(5  30  +  ^  C$cQ  (  S$Ct(H-C$S0S^  ) 

w  ■ 

“y 

- 

9c9  -  $isdc6  (  c$ct|/-s$s8si|> ) 

s  $  +  <JJS0 

It  is  well  known  that  for  some  applications,  for  example,  for  the  tethered  Shuttle 
subsatellite  system,  the  mass  of  the  Shuttle  is  much  greater  than  that  of  the 
tethered  subsatellite,  so  the  center  of  mass  of  the  Shuttle  can  be  considered  to  be 
the  mass  center  of  the  whole  system  and  the  shift  of  the  center  of  mass  of  the 
system  can  be  neglected  i.e.  a„-f#=«0  in  equation  (6).  However,  in  our  system  the 
shift  of  the  center  of  mass  of  the  system  will  be  considered  and,  in  general, 
a«-f<>¥0;  it  will  be  calculated  in  the  development. 

After  projecting  equation  (6)  on  the  tether  frame.  Rt.  the  following  is  obtained 

••  «  • 

(ag-fg)  |  c+TCr+2ojxr-Huxr-Kjirx(uxr) -Mr]da-  T-(?)[c-«]cdm-0  (8) 

where  [,  indicates  the  projection  onto  the  frame.  R.. 

The  expression  of  the  r  of  the  tether  system  is  different  from  that  of  the  shell 
due  to  relative  motion  of  the  tether,  so  we  consider  the  tethered  subsatellite  system 
and  the  shell,  separately. 

Tethered  subsatellite  system 

r  =  T'1  r,  +  h  -1-  q,  (9) 
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where  7*  -  (  x  +  u,  v,  w  )  (10) 

is  the  position  vector  of  dm  from  O  projected  onto  the  frame,  R„  u  represents  the 
tether's  longitudinal  elastic  displacement,  v,  vv  represent  the  displacements  in  the 
orthogonal  directions  transverse  to  the  OX  axis. 

h  =  (  h„  hv,  ht)  is  the  position  vector  of  point  O  from  the  shell's  center  of 
mass,  O,. 

q„=»(u,#(  0,  0)  is  the  shell's  elastic  displacement  vector  of  the  apex  of  the  shell 
(according  to  the  assumptions  there  is  only  elastic  displacement  along  the  X,  axis). 


Hence, 


r  -  (T'1)  r,  +  qa  +  (T")  r. 


r  =  (T-1)  rt  +  2(T*1)  r,  +  (T")  r.  -  q, 


0  mu  oj 

z  y 

0)  0  mu 


mu  u 

y  x 


According  to  vector  algebra 


where 


'jjxr*ux(<uxr)  -Mr  *  [QJ  r 


Cq]  *  C«]  +  Du]Du]-£m] 


After  substitution  of  equations  (11H14)  into  equation  (S)  there  results 

[fa,  -  f,  )  ,  -  a.j  dm  -  L(q)  -  el.  dm  tin) 

where  a,  -  r.  +  Tq„  +  2[TfT"1)  -1-  T[  a  ]T'1]  rt  -  2T[  y  ]q,  + 

[T(T"1)  -  2T[  y  jff")  -  T[Q]T'']  %  m  TO}(h-q,)  <lf>) 

In  order  to  form  the  svstsm  linear  cuuaiion<  of  motion  cu nation  (lf>)  is 


linearized  as  foilows: 
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(17) 


at  *  (at  ,ae  ’ac  5 
x  y  2 

a  /&r-4"+M"-u  "-3u+3u  -2w'+<2a,+28,-3)x 

x  c  p0  p0  _ 

+■  0-20  ')h;t+<({)',+2^'-4<l)+Y)hy+<30-e,,)h2 

ac  /aJc  "  v''4v^<Y,,-WYH>,,-4(ii)x+(c|»',+4<$i-3Y)h  +h  -(^"+^)h 
y  x  y  z 

a  /coJ-w"+2(2,-Hi,)-2u:  -(a"+3"+30+3a)x 
z  c  p0 

+<8"+3a+30)h  -K2$'H/"-H0h  -28 'h 
x  y  2 

where  (  )'  *  d(  )/dx  and  t  =»  t,  2  is  the  length  of  the  tether 


Shell  system 

Now  consider  equation  (8)  for  the  shell  system 

r  =  7,  +  q,  =  (  x,  +  u„  y„  z,  )  (19) 

7  -  (u„  0,  0)  ;  7  -(ii„  0,  0)  (20) 

where  x„  y„  z.  are  the  coordinates  of  dm  of  the  shell  in  the  frame,  R„ 

After  substitution  of  equation  (20)  into  equation  (8)  there  results 

[(a,  -  f„)|t  +  a,  ]  dm  -  T  L(q,)  -  e  |(  dm  =  0  (21) 

where  a.  =  [3u,  -  u„  0,  -2u,']T  +  T  [Q][x..  v..  z.  ]r  (22) 

Equation  (15)  is  integrated  over  the  tethered  subsatellite  and  equation  (21) 
is  integrated  over  the  shell.  The  two  results  arc  added  together  and  it  is  obtained 


(a,  -  f,)1,  =  (i  m_)  [  E  -  /  a.  dm  -  /  a,  dm  ]  (23 

s.c  ? 

where  mT  -  m,  m*  +  m,  :  E  »  /  7  dm 

s. 

and  m,.  m..  m,  are  the  mass  of  the  subsatcillrc  .  ■  c mb  shell,  respectively.  Tin. 


subscript  s.t  or  p  designates  the  integration  over  me  Nid-m.-ucilite  anti  tether  or  the 
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shell,  respectively. 


To  obtain  the  Rayleigh-Ritz  solution,  u,  v.  w  can  be  expanded  in  series  form 

in  terms  of  a  set  of  admissible  functions1*' 

u  -  Z  (x)A  (c) ;  v  -  E  D  (x)B  (c) 
n  n  a  ti  ti 


w  »  E  <t>  (x)C  (C);  L(q)  -  [L(u)  ,L(v)  ,L(w)  jT 
n  a  a 


Introduce 


I,  »  '  $  dm  :  I.  »■'  'U  da  :  I  *  ^  xdm 

*  s.c  a  ’  3 . c  ’  x  3,c  050 

a  a 

V.  *  •'«  V.-  V.  ’si  Vn1* 

8»  •  ,{e  x2dB  :  Vf.fe  ’V-  :  *»„»,*  dcW" 


Hence 


/  udm  «  El,  A  ;  J  u'dffl  »  El.  A' 

.  _  a  a  ~  ’S'  a 

3 ,  c  a  s.c  a  a 

/  u"dm  *  El.  A"  ;  /  v  da  »  E  I.  3  ,  ecc  nt\ 

a  '!>  m  a  $  a  Uo) 

s.c  a  s.c  a 

Mathematical  expressions  for  the  natural  frequencies  and  mode  shapes  of  the 
transverse  vibrations  of  a  shallow  spherical  shell  with  a  completely  free  edge  have 
been  obtained  in  Rcf.[9]. 

u  *  EA  (c)4>  rpT) 

P  n  Pn  ?n  (Z/J 

where  $  is  the  nth  mode  shape  function 

*pn  "  ^ j «s7 K,+,J"x<' ^ j  k’  ^  +  Djk:Ik(Ajk.’^cosk(34^0)  (28) 

a  is  the  base  radius  of  the  shell 


Since 


f  2  da  *  0 
P  ?ti 


hence 


*  u  da 

4  3 

3 


;  u '  da  *  >  u"  da  *  0 

3  3  3  3 


[t  is  expected  that  the  natural  'Ycauencic^  and  'h noe  •'unctions  will  be  modified 


bv  the  presence  of  the  tether  s  vs  tern.  However,  m  ->u  nor  wc  still  adopt.  3. 


as  the  assumed  shape  functions,  since  the  tether  system  mass  is  much  less  than 
that  of  the  shell. 

Equations  (23)  ,  (18)  ,  (22)  ,  (26)  ,  (29)  arc  substituted  into  equation  (15) 
and  the  resulting  vector  equation  is  projected  along  the  X,  Y,  Z  axes.  After 
integrating  the  projection  along  the  X  axis  over  the  subsatellite  and  the  tether  the 


translational  equation  for  the  tether  motion  can  be  obtained 

-  -  cr'-u”  +3u  )+EI* \a"-3AJ  -  21^*  C '+<2a  ’+29  ' -3)  I* 

cx  c  sc  pn  p-  m  il  a  a  a  d  a  '  x 

u  u  a  a 

+»*eC(3-2e>)hx+  (*"+»  Bstx*  sp  s; 

where  *sc  -  ;  n*c  -  aje*p/«s  ;  I*  -  a?y»s 

I.t  *  3  I .  /m.  ;  I*  •  a  I  /a. 

"a  p  "a  “  Ja  _  p  *a  ‘ 

E,„,  E„  are  components  of  E,„  E.  along  the  X  axis 

E„  =»  /  e  dm  ;  E ,  -  f  c  dm 
s.  c  P 

and  F.„  is  the  tether  tension 


The  equations  for  the  rotational  motion  of  the  'cthcrcd  subsatellite  can  be 

obtained  by  the  following  operation 

/  r,  x  equation  (15)  =0  (33) 

s.  c 

By  projecting  equation  (33)  along  the  Y  and  Z  axes,  respectively,  the  rotational 
equations  for  the  pitch(  in-plane  swing  )  and  mil  i  oin -m'-pianc  swing  )  motions  are 


obtained  as 


o"  -3  )  -  3Z(h*  -a  :*  )3  -=• 

a  a  a  a  :o  ::  a  a  :c: ' 

a  a  a 

(.i  -a  -  )(j  -«o-3v)  a  ;  a  .J  —  v.  ... 

:cc  x  :c  x  a  a  a  7  a. 


.  -  a 

ia  c 


(34) 

(35) 
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where 


lA  -  H  .  -I  I  /m. 
x$_  x<b  x  4  Z 
a  a  a 

*x*  *  a*  -V*  /nz 

a  a  a 


H 


xx 


H 

XX  X  I 


(36) 


Ltv,  L„  are  components  of  the  torque  L„  produced  by  the  external  force 
By  the  following  operations, 

{  'K  [Eq.  (15)  ]_»  0  (37n); 

(39  n) 


s.c  a 


,{C  (,S'>  3/  0 


(38n) 


/  *  CEq.  (15)3-0 
s.c  a  2 


the  nth  longitudinal  and  vibrational  mode  equations  ar?  obtained  as 
I*  +3u  )+  t  a*  „  (A" -3 A  3-201*  C'+(2a'+20  '-3)1! 


i>0  “  V."' 


J  "  I  i  w 

a  m  v  v  a 
a  a 


xi|» 


+r  c ( 3-20  ')h  «W+2$’-4<fr+Y)h  +(9,,-30)h  ]+  z:<  A  -  H(n) 


a  an  a  ex 


31  a  m  n 

I*  C(4.r'+4*-3Y)hx+hy-<*"-H(,)h2]+H^^Uba  ■ 

EH*  .  C"+2I*  (1 ' -u ’  )+2EH*  .  A'  -H*  (a"+e"+3ct+3e) 

*  4^  a  in  pQ  y,  a  x<fia 

+1*  C(eM+3a+30)h  +<2p’-H4M>")h  -29 'h J+H.  .  -Jc  =■  H 
i  x  /  2  +>  ts  n  a 


where 


a  a 


(n) 

ez 


n  a 


am  a  a 


a  a 


a  a  a  a 


(40n) 


(4  In) 


(42n) 


(43) 


*  m*pt  -ri  S  /:aE 
am  nym  a  a 


and  the  terms  H„  .  Hfy,  H„  results  from  the  external  force.  Z  -^A  ,  '.*£3  , 


"> 

arc 

a  a 


result  from  the  elastic  force1 ! 


The  equations  of  shell  motion 
By  the  following  operation 

/  r„  x  equation  (fi)  =  0 

p 

the  equations  of  the  rotational  motion  of  'he 
ea  nation  (M  includes  the  tether  force  acting 


^ It e ■  •  m  obtained  as  (note:  c  in 
on  .no  " no: I  which  can  be  obtained 
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by  integrating  equation  ( I V  and  it  represents  the  effect  of  the  motion  of  the 
tethered  subsatellite  system  on  the  shell). 

-  a/J*){ffl*chyh2  /*T  -hy[l*cc“+(I*-hx3*sc)(9"+38) 


-2a*  ( i’-u’  -h  8’) -211*  A’ -SI*  G“]+h  Cl*(Y"+Y) 
sc  Pn  z  i|<  a  a  nr*  2U  xv ' 

,  u  —  a  a 


-<I*-h  a*  )  (<f>"+4<$)+EI*  (3"+B  )]+<L_  +  L  )  /cu2  } 

x  x  sc  a  $  a  a  tpx  epx  c 

a 

*-(i/jp(3b1c>^vt5/.*+biC»;t(»-,jo*3«Po) 

ta  ‘  ni  ^ 


(44) 


)+2II-I  A,«-I*a"+hya,«(»’+*-|’")><LSp7+l.p^/"c  >  (45) 

U  a 


•  C-1/-T*)(h  h  -3I*)/ca2*h  .CtY,,+v):*-a*  (V’-HOh 

2  *•  5  /  SC  X  «  A  X  * 

+SI*  (3"+B  ) ]-h  [a*  (Z”-u”  )  *?( I*a’+<l*-m*  h  )8'-SI*  C‘) 

a  «  ®  y  sc  pn  ?n  x  x  sc  x'  m  4  a 


*4<39-S")VEI,;  )/^> 

m 

where  S*  -  /J*  =  *  <VV/Jy  ;  ns  * 

J*  *  J  +a*  (h2+h2)  ;  J*  *  J  +m*  (h2+h2)-I*'n 
x  x  sc  y  z  *  y  y  sc  z  x  x  x 


(46) 


J*  -  J  +m*  (h2+h2)-I*h  ;  I("y 
z  2  sc  x  y  x  x  I 


/  x  p  dm  ;  C  »  A_  /a 
p  ?  Pn  ?a  ?a 


(47) 


L«„,  Lepv,  L€oi  are  the  components  of  torque,  produced  by  E,t  and  E„  which 


appear  in  the  tether  force  acting  on  the  shell:  L ......  L.ov.  L„t  are  components  of 

torque  which  are  contributed  by  the  external  forces  -'cong  or  the  shell:  and  Jx,  Jv, 
Jt  are  the  principal  moments  of  inertia  of  the  undeformed  shell. 

Bv  the  following  ODeration 


/  o  [  equation  (6)  ],  -  f> 
p  ’  - 

he  nth  shell  elastic  vibrational  mode  ecnation  ;s  om-med 


-v-;-3)e  *2iy‘h’/:^3.  -  r*'- ?  . 


o 

■  n 


-E  ,-a  „  / M  a 
x  -d  a  :  n 


(48) 


(49n) 
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where  E,  is  the  modal  component  of  the-  external  force.  are  $•-  at  the  point,Of 

?a  • 

M,  is  the  nth  modal  mass,  F„  is  the  component  of  the  tether  force  acting  on  the 
shell  along  the  X.  axis, 2  /ci»„  where  cuu is  the  naturnl  frequency  of  the  nth  mode. 

n  w  rTl 

Equations  (30)  .  (34)  ,  (35)  .  (40n)  .  (4 In)  ,  (42n)  ,  (44)  ,  (45)  ,  and  (46), 
(49n)  compose  the  complete  system  equations  of  motion. 

Now,  for  our  special  case  (Fig.3).0  is  along  the  shell  yaw  (i.e.,  X,)  axis,  hence, 
hy^h^O  and  it  is  assumed  that  there  are  no  external  forces  acting  on  the  system. 
By  examination  of  the  equations  for  this  special  case  the  following  conclusions  can 
be  reached:  (1)  the  shell  roll,  yaw  motion,  tether  our-of-plane  swing  motion  and 
elastic  vibrations  are  decoupled  from  the  shell  pitch  motion,  shell  elastic  vibration, 
tether  in-plane  swing  motion  and  in-plane  elastic  vibrations:  (2)  the  shell  pitch  and 
elastic  motions  are  coupled  directly  to  each  other  through  their  rates;  (3)  since 
I,"”  =»0  for  ail  shell  elastic  modes  except  for  the  axisym metric  modes,  only  the 
axisymmetric  modes  are  coupled  to  the  shell  pitch  motion:  nonaxisymmetric  modes 
are  independent  of  the  system  motions,  and  would  have  ro  be  controlled  separately 
within  the  linear  range. 


2.3  Stability  Analysis 

It  is  well  known  that  the  shell  pitch  and  roll-yaw  morions  are  unstable  about 
the  present  nominal  orientation  as  J,  .J.'-J,  without  r he  attached  tether  system. 
In  the  present  paper  stability  conditions  for  the  ’chercd  shell  system  will  be 


deveionea  wnen 


:s  considered 


tonsidcred 


riaid). 


In  generni.  a  finite  number  of  elastic  modes  in 


he  mooei  :s  -.o  he  retained  for 
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practical  purposes  (truncated  model).  In  the  present  paper  a  few  such  truncated 
modes  are  considered. 


Rigid,  constant  length  tether  for  in-plane  motion 

In  this  case  all  of  the  tether  elastic  modes  arc  neglected  and  the  tether  length 
is  fixed  (without  tension  control).  Hence,  according  to  equations  (34),  (45)  for 


our  special  case  the  equations  of  in-plane  motion  arc  simplified  as  follows: 


X1a,,+6"+3a+38 


:<-a',+9”-3ft  3 
2  y 


where 


*  .  ,  *  it  ,  * 

X,  -  H  / (H  -I  h  );  K  -  -I  h  /J 

1  XX  XX  X.  X  2  XX} 

a*  -  (J  -J  -«*  h2s-I*h  )/J* 

y  x  z  sc  x  x  x  y 


The  system  characteristic  equation  is  given  by 
(K1-K2)X4+3(l-iCLa*-X2)X2-9n*  -  0 


*.  * 

y 


since  K  -X-  -  Ch  J  +h  a  a  (a  +a  /4)/3m^]/  (H  -h  I  )J 

xx  y  x  c  ?  s  t  -  xx  x  x  j 

if  h„  <  0  then  K,  -  K*  >  0 
The  neutral  stability  conditions  are 


n*  <  o  (55):  i >  o 


y  %2 


Meanwhile,  if  h,  >  0  the  neutral  stability  conditions  ire 


K,  -  K-  >  0  and  r.'  <  0 


'■'ut  (58)  is  inmost  I m possible  to  satisfy  if  %  ^ 


(50) 

(51) 

(52) 

(53) 

(54) 


(56) 

(57) 


and  9(1-X,.T  -K„)"  +  4(X,-X-)9Q  >  0 

-y  ii/ 

It  can  be  proved  that  if  condition  (55)  is  satisfied,  then  (56),(57)  arc  also 
satisfied. 


(58) 


'ooordinu  to  (52>  and  •  >5> 


:t  is  setter  to  cnoose  a. 


0.  so  that  the  aeu'ra:  o'"n\  conditions  for  in-niane 
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motion  are 


h*  <  0  or  h^-h,  >  0 

J,-J,  +  m,tTi2  +  I/h«Jt-J,+[(m,  +  mt)hs  +  (mt-)-m..2)h  ]  m,/m  >0  (59) 


Rigid,  constant  length  tether  for  out-of-plane  motion 

According  to  equations  (35)  ,  (46)  and  (44)  the  equations  for  the  rigid, 

constant  length  tether  for  out-of-plane  motion  are  simplified  as  follows: 

rC3Y',H>H+K3Y-4i2*c!)-(l-a «  0 


'r-flV(w*)4'  -  o 

X  X 


where 


X 


-I*h  /J* 

XX  Z 


3*  »  J  -J  • +(a*  h2-I*h  )/J* 
Z  >  X  SC  x  X.  xy  z 


n 


X 


where 


l0  -  :VK3 


The  system  characteristic  equation  is  given  by 

6  4  2 

a„A  -Ha-A  +a,A  *  0 
0  2*6 

-K, 

a2  “  “flx(Sl"K3J+3+iCi*WiClas  ) 

a4  “  ^\(3^^*’^-3K3)+4r**(3+K1) -4K^-(3+K1)(l-np(I+nx) 

a6  *  -^xC(3-«i)n*-K3] 

the  neutral  stability  condition  for  this  system  is 

a2  <  o 

i.e.,  \-  is  a  negative  real  number. 

Since  a,  >  0.  condition  (64)  is  counted  !n  rho  ni'nwlna  conditions 
4„>0  (65)  ;  a-,a.-a.a,>0  (fi6)  ;  a, >0  (^7) 

-  4.  ■*  J  0  J 


(60) 


(6t) 


(62) 


(63) 


(64) 
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and 


a  -  a^aj4.iaa0a2a4a6-4a0aj-27aja^-4a6a^>0  (68) 

The  conditions  (65)-(67)  a  re  the  Routh-Hunvitz  conditions  resulting  from  the 

cubic  equation  in  the  variable  X2  [equation  (62)  ]  and  condition  (68)  is  from  the 
condition  that  this  cubic  equation  has  three  real  roots. 

It  is  known121  that  one  of  the  necessary  and  sufficient  stability  conditions  for 
out-of-plane  motion  of  the  shell  system  without  the  attached  tether  system  is 

ax  "  (J2“Jy^/Jx<  0  <69) 

Now  we  assume  condition  (69)  is  still  satisfied,  so  from  condition  (67)  it  is 

obtained  that 

Q*  >  :<3/(3+K1)  (70) 

If  the  condition  (70)  is  satisfied,  it  can  be  proven  that  the  conditions  (65)  and 
(66)  are  also  satisfied  and  it  is  demonstrated  numerically  that  condition  (68)  is 
also  satisfied  for  a  variety  of  system  parameters. 

Hence,  the  neutral  stability  conditions  for  out-of-pianc  motion  arc 

\<  0  (71)  ;  r-2f>:<3/  (K1-3)  (72) 

The  neutral  stability  conditions  for  a  rigid,  constam  length  tether  are  conditions 
(71)  .  and  (72)  together  with  condition  (59) 

Typical  stability  regions  for  in-planc  and  out-of-plane  motion  in  the  parameter 
space  m„  h  are  shown  in  Fig. 5:  it  is  seen  that  die  suability  region  for  in-plane 
motion  is  larger  than  that  for  the  out-of-pianc  motion. 

It  can  be  proven  chat  the  out-of-planc  motion  win  asymptotically  stable  when 


lamotna  or 


tin  tz:c 


loiiei  nc- 


Imons 


r:.) 
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Rigid  variable  length  tether  with  Rupp's  tension  control  law  for  in-plane  motion 
According  to  equations  (30)  .  (34)  .  (45)  the  equations  of  motion  with  tension 


control  are  simplified  as 

e”+2Ka,+<K+3  )28'-3Ks-  F  /w*  »*  l  +3  (K+S  )  -  if 

X  CX  C  St  * 

X1a,,-H”+3a+3e-K4e'  -  0 

K^a  ''+8  "-38*9  -iC^e  '  -  0 

where  3  »  h /l;  e  »  (i/z  )-i  ; 

•'  c 

K,-  21*4/ (H*  +hl*)  5 

4  X  xx  X 

For  Rupp's  tension  control  law1 101 

Af  -  -<K  e+K  * s  * )  (75) 

£  £ 

the  system  characteristic  equation  is  developed  ami  the  Routh-Hurwitz  criteria 
applied.  After  some  complicated  algebraic  manipulations  the  following  expressions 
for  the  principal  minors  are  obtained 


(73) 

K*  (aig-Hm^/2)  /  (ras+m  ) 

K  -  2m*  h£/j*  (74) 

5  su  y 


°1  *  r:;i'K2)Kr' 

°2  •<:<1-K2)Ke .  Cc^-k^k^xcx^-^k^-Xj)  +2Sx(K5X1-K4:<2)] 
D3*(:<1-K2)xJ,C-3^(:<3:<1-K4:<2)(2:<1-s:;4)3:c+6K(x3:<1-x,:<2ric4-x3)2/x.] 

D,  »9(:<1-x2)k^ , { (A*a*2+a*a*+c*)x*  kskci-k^/Xj)  C(X4:<2-X5:<1)^ 

-<VV  (sl5k1-k4k2«4-ks)  3  2> 


i* 

y 


D-  *  -64ar.*x^x(i-x1)C(x4x2-x5:<1)n7  -x5(x5:<1-x4x2^c4-x3)/x4]z/x5 


-9n*x*D. 

y  s  a 


(76) 


where 


“Jrv 


ta,' 
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8*  -  8Kh2 4  V,»t (  VV4>  [3or CB^+tll*) *J*2  ] 

C'  *  '■I*«CJ!r4m;i.h2)j2/C(H^+hI*)J*3l  (77) 

The  necessary  and  sufficient  stability  conditions  arc 

D,  >  0  (  i-l,  .  6  )  (78) 

Since  h  >  0  (  or  h.  <  0  )  according  to  equations  (54)  and  (51) 

K,  -  K,  >  0,  K,  <  1 

Hence,  from  Di  >  0,  D,  >  0,  and  Ds  >  0.  there  results. 

X  .>  0,  s*  >  0  and  a*  <  0  (79) 

s  z  y 

Meanwhile,  since 


and 

thus,  if 
then 


’  21^ZJy/^(a«+hIx)JP  >  0 

:<r<l-K4:<2  »  2hi.3a  a  (a  -ha  /4)/  [3(H*  +hl*)  J*ju]  >  0 

r  ^  U  AA  X  j 

2Xx-KX4  -  2A2m6(«#+Bic/4)/C3(aJsc+hI*)msC]  >  0 
3*2  -  4A*C*<  0 

x_,>  o,  x!>  o,  n*<  o, 


D(  >  0  (  i-  l. 


6) 


(80) 


So  the  necessarv  and  sufficient  conditions  for  in-nlanc  motion  stability  with 


Rupp's  tension  control  law  are: 

2*  <  0;  X. ,>  0  and  X*  >  0 


(81) 


Stability  conditions  for  flexible  tether 

When  tether  flexibility  is  considered  it  is  difficult  to  qer  analytical  results  for  the 
stability  conditions  because  of  the  hish  order  of  the  s-.ccti.  However,  it  is  found 


numerically  that  the  system  stability  conditions  are  Mechanized  both  for  in-nlane 
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motion  (  with  constant  length  tether  or  with  Rupp's  tension  control  law  )  and  out- 
of-plane  motion  (  with  constant  length  tether  )  when  the  tether  flexible  modes  are 
included  in  the  system  model,  for  a  variety  of  system  parameters  (  for  example,  for 
variations  of  h,  m„  and  l  ). 


2.4  Optimal  Tension  Control  Law  for  In-Plane  Motion  during  Stationkeeping 

\ 

In  the  last  section  it  has  been  demonstrated  that  the  system  is  asymptotically 
stable  with  Rupp's  tension  control  law  for  in-piane  motion.  However,  in  order  to 
improve  the  transient  responses  two  alternate  optimal  control  laws  are  introduced. 
One  is  an  optimal  control  law  based  on  tether  length,  in-plane  swing  angle  and 
shell  pitch  angle  and  their  rates  for  the  rigid  massive  tether  model;  another  is  an 
optimal  control  law,  which  includes  additional  feedback  of  the  tether  vibrational 
modes  and  their  rates.  For  the  system  with  the  state  variable  format  equations  the 
optimal  control-  U,  which  minimizes  the  performance  index 


J  -  /  (  XTQX  +  UTRU  )  dt 

o 

is  given  by  U  -  -  (  R“BrP  )  =  -KX 
where  X  is  the  state  variable 

Q  is  a  positive  semi-definite  state  penalty  matrix 
R  is  a  positive  definite  control  penalty  mntrx  and 

P  is  the  positive  definite  solution  to  the  steady  state  Riccati  matrix  equation 
-PA  -  A'P  -  P3R-'3rP  -  Q  =  <*• 

The  ootimai  control  iaw.  based  'in  the  'ether  ienirn.  •.'.-plane  swing  angle,  and 
sneil  pitch  ansie  and  their  rates  for  the  riaid  tether  ••mnici  :  with  state  variables 
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XT-  [ 


a,9,e,a’ ,3 ' ,e 


])  takes  the  form: 


<lf  -  -<Kaa'+Kg3+S££+Kci,ci,+K3I9,-HCr  ,e’) 

For  the  system  parameters:  m,  -  10000kg,  m,  -  500kg.  -I  km,  m,»  8.35kg, 
h^O.OSkm  some  typical  system  simulation  (  three  vibrational  modes  are  included) 
results  show  that  the  transient  responses  for  the  rigid  tether  optimal  control  law  is 
better  than  that  for  Rupp's  tension  control  law  for  some  of  the  control  gains, 
especially  for  the  damping  of  the  tether  and  shell  pitch  angles.  But  it  is  also  found 
that  the  system  could  be  unstable  for  some  of  control  gains  (  Table  1). 

It  is  obvious  that  the  transient  responses,  when  the  optimal  control  law  includes 
feedback  of  the  vibrational  modes,  are  better  than  rhe  responses  based  on  the 
previous  control  law.  The  improvement  is  especially  noted  in  the  damping  of  the 
tether  vibrational  modes.  A  typical  comparison  of  transient  responses  for  the  three 
different  tension  control  laws  is  shown  in  Figs.  6a.  6b.  and  6c. 


2.5  Conclusions 

The  orbiting  shallow  spherical  shell  pitch  and  roll-yaw  motion  are  unstable  when 
the  symmetry  axis  nominally  follows  the  local  vertical.  However,  it  is  suggested  that 
gravitational  stabilization  could  be  achieved  by  including  a  tethered  subsatellite 
system  to  provide  the  favorable  moment  of  inertia  distribution.  The  tether  could  be 
connected  at  the  end  of  a  rigid  boom  which  is  attached  to  the  shell's  apex.  The 
equations  of  morion  for  such  a  tethered  shallow  spi'c'-ctd  shell  in  orbit  with  the 


-i resent  nominal  orientation  na*  c  eec.n  oeveiorc::  m 


■'■er. 


i  ne  sneii  roii-vaw  motion,  tenter  out-u-rinne  >wm<z  moron.  and  the  tenner  out- 


of-plane  elastic  vibrations  are  decoupled  from  the  shell  pitch,  shell  elastic  vibration, 
tether  in-plane  swing  motion  and  tether  in-plane  elastic  vibrations.  For  given  shell 
and  tethered  subsateilite  system  parameters  a  suitable  rigid  boom  length  could  be 
chosen  in  order  to  provide  a  gravitational  stable  structure  both  for  in-plane  and 
out-of-piane  motion.  The  in-plane  motion  of  the  system  could  be  asymptotically 
stable  with  Rupp's  tension  control  law.  It  is  demonstrated  numerically  that  the 
flexibility  of  the  tether  would  not  affect  the  stability  conditions  for  the  constant 
length  tether  or  for  the  variable  length  tether  with  Rupp's  tension  control  law  for 
a  variety  of  system  parameters.  The  transient  responses  can  be  improved 
significantly,  especially  for  the  damping  of  the  tether  and  shell  pitch  motion,  by  an 
optimal  control  law  for  the  variable  length  tether  model.  It  is  also  seen  that  the 
system  could  be  unstable  when  the  effect  of  tether  flexibility  is  included  if  the 
control  gains  are  not  chosen  carefully.  The  transient  responses  can  be  further 
improved  by  including  the  state  feedback  of  the  tether  vibrational  modes  into  the 
optimal  control  law,  especially  for  the  damping  of  the  tether  vibrations. 

Extensions  to  the  present  paper  could  consider  the  effect  of  the  shell  flexibility 
on  the  system  stability  and  control  and  some  kind  of  active  control  could  be 
introduced  (in  addition  to  tether  tension  control!  to  improve  system  performance. 
Additional  control  will  be  required  to  provide  for  out-of-plane  damping  of  rigid 
motions  and  vibration  suppression. 


Table  1.  Stability  Characteristics  for  Different  Control  Gains 


R 

K 

a 

H 

Kc 

Xa-' 

V 

V 

Stability 

1 

6.938 

5.875 

6.114 

2.722 

m 

4.917 

unstable 

2 

5.456 

4.841 

6.034 

1.637 

2.455 

4.377 

stable 

5 

4.296 

4.059 

5.984 

0.583 

1 .058 

H 

stable 

10 

3.812 

3.742 

5.967 

0.012 

0.304 

3.478 

stable 
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3.536 

3.565 

5.959 

-0.405 

-0.249 

3.217 

stable 

30 

3.437 

3.503 

5.956 

-0.591 

-0.499 

3.095 

unstable 

where  Q  - 

1  Jij 

Table  2. 

Control  Gains 
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4.84 
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0  0 
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Fig.  3.  Tethered  Antenna; Reflector  System 
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Fig.  4.  Coordinate  Systems 
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h  >  h  stable 


Fig.  5.  Stability  Regions  for  fn-Plane  and  Out-oi-Pinnc  Motion 


a  (deg) 


9  (deg) 


Fig.  6a.  Transient  Responses  for  Rupp's  Control  l.aw 
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Tiaie  (orbics) 

Fig.  6c.  Transient  Responses  for  Flexible  Tether  Model  Optimal  Control  Law 
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3 


REVIEW  OF  THE  CONTROL  OF  TETHERED  SATELLITE  SYSTEMS 


This  chapter  reviews  the  steps  ir.  the  development  of 
control  laws  for  Shuttle/Platform-Tethered  Subsatellite 
systems.  The  tethered  subsatellite  systems  have  been 
proposed  for  numerous  applications .  This  has  led  to  many 
investigations  dealing  with  the  dynamics  and  control  of  such 
systems  during  their  deployment,  station  keeping,  and 
retrieval.  A  brief  comparison  of  control  laws  used  by 
various  investigators  is  described  here  in  order  to  evaluate 
different  control  methods  for  tethered  subsatellite  systems. 

Finally,  recommendations  are  made  as  to  the  suitability 
of  the  different  control  laws  for  adaptation  with  the 
proposed  orbiting  tethered  reflector  systems. 

3.1  Deployment 

This  operation  involves  moving  the  subsateliite  from 
the  Shuttle / Platform  to  a  distance  as  much  as  100km  from  the 
Shut  ;ie/Platf orm .  Since  most  of  the  useful  mission 
activities  start  after  deployment  of  the  subsateliite,  it  is 
desirable  to  deploy  the  subsateliite  to  the  cperationai 
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mass,  ms ,  attached  to  one  end  of  the  tether  is  ejected  from 
the  Shuttle  and  performs  a  free  flight  until  the  tether 
becomes  taut  at  which  time  the  payload  is  subjected  to 
impulse  affecting  the  motion.  Impacts  and  free  flight  occur 
alternately  until  enougn  energy  has  been  dissipated  during 
impacts  so  that  no  further  ones  occur.  The  tether  then 
becomes  permanently  taut  and  the  system  behaves  like  a 
spherical  pendulum.  If  seme  viscous  damping  is  provided, 
the  pendulum  settles  along  the  local  vertical.  Whether 
deployment  in  the  desired  direction  (upward  or  downward) 
c  :curs  or  net,  the  duration  of  deployment  depends  on  the 


initial  ejection  velocity. 

Deployment  is  mere  likely  to  be  carried  out  using  an 
active  control  system  to  guarantee  adequate  dynamic 
performance .  Rupp  [10]  first  made  a  preliminary  treatment  of 
“he  dynamics  of  the  Shuttle-Tethered-Subsatellite  system  in 
which  the  motion  is  assumed  to  occur  only  within  the  orbital 
plane  and  tether  mass  is  neglected,  and  set  up  a  tension 
control  law  m  the  form 

T  =  KiL  -  1 1 L  *  K:lc  (82) 

where  L  and  L  are  instantaneous  length  and  length  rate, 
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•c:\trol  law  requires  the  measurement  of  tether  tension,  rate 
f  ohe  r  d*  cl  fnter.t  and  instantaneous  tether  length. 
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tional  contact  with  the  tether.  The  tension  can  be  measured 
by  a  spring  damper  arrangement  on  the  same  pulley.  The 
measurements  are  fed  into  a  computer  which  calculates  the 
required  torque  that  must  be  produced  by  the  motor  driving 
the  tether  reel  system. 

In  the  dynamic  simulation  using  Eg.  i'll;  Rupp  [10]  used 
K2  =  K]_-3ms^-  while  two  values  of  K;j_(5msto-  and  7msw-)  were 
chosen,  [co  is  the  assumed  (circular)  orbital  frequency  of 
the  Shuttle  orbit.]  The  coefficient,  Cj,  was  used  to 
critically  damp  the  longitudinal  stretching  oscillations . 

The  control  law  was  quite  effective  in  damping  in-plane 
motion  during  deployment,  but,  for  out-of-plane  motion  was 
not  considered. 

Eaker  et  al .  [12]  treated  the  three  dimensional  dynamics 
and  control  including  the  inertial  effect  of  the  tether 
mass,  aerodynamic  heating,  and  aerodynamic  forces  on  the 
tether  and  sub satellite.  The  form  of  the  control  law  used 
was  the  same  as  in  Rupp  [10]  with  the  exception  of  modifi¬ 
cation  cf  the  commanded  length  for  desired  deployment  and 
re crieval  maneuvers.  The  modified  tension  la w  was  in  the 
form  as  follows: 


m  + 
s 


:  r.-~c  ;  co-i. 


!  83 


where  m  .•  and  ms  are  the  mass  cf  subsatel 
respectively .  R  is  the  ratio  between  the 
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damping  ratio.  A  commanded  length  of  the  form 

Lc  =  KXL  +  K3  (  84) 

was  suggested. 


The  combinations  of  exponential  and  uniform  rates  of 
change  of  length  can  be  expressed  in  the  form 


'Li  exp  (ct) 

Li< 

Lc  <  Li 

r*  * 

Li  (lect) 

VI 

Lc  <  l2 

85) 

L3  -  L3 -Li  exp  (-ct) 

V  1 

c  i 

•n 

Lc  1  Le 

where  c  is  a  positive  constant,  and 

Lf  are 

initial 

and 

final  lengths,  respectively,  while  L\  and  L3  are  two 
intermediate  lengths. 

Deployment  is  basically  a  stable  operation  [Baker  et 
5.1 .  )  L  l'2i  ;  however,  towards  the  end  of  deployment  ,  when 
aerodynamic  effects  become  important  bounded  steady  inplane 
rotational  and  elastic  oscillation  may  result.  In  addition. 
•:  ut-c  f -plane  rotaticns  and  vibrations  recur  for  eccentric 
oroits  inclined  to  the  equatorial  plane.  This  motion,  does 
net  cav.se  any  serious  problems  and  can  be  eliminated  once 
the  deployment  is  ccmp-eted. 


Kaiaghan 


et  ai .  [13]  used  a  finite  element  approximation 
ec  mathematical  model  of  the  Shuttle-Tethered- 
ystem  inch  inc_udes  tether  mass  effects; 
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L '  =  Lc  '  L 1  J  (86) 

where  K  and  X  are  gains  and  state  vectors,  respectively. 

Subsequently ,  Bainum  and  Kumar  [  15J  first  developed  an 

optimal  control  law  based  on  an  application  of  the  linear 

regulator  problem,  with  control  provided  only  by  modulating 

the  tension  level  as  a  function  of  the  difference  between 

the  actual  and  commanded  length  (e  =  Lc/L  -1),  actual  length 

\ 

rate  (£ ' ) ,  inplane  tether  line  swing  angle,  a,  and  its  rate 


a  '  ,  i  .  e , 

T  =  K-e  +  K  £  '  +  K  a'  -c  K  a  '  -r  Tr  ;  T0  =  3a?L„ms  (37) 
c  £  a  a*  c  '  o  -s'-*'' 

The  system  was  idealized  as  two  point  masses  connected  by  a 


massless,  mextensible  tether  with  the  system  moving  in  a 
nearly-circular  orbit.  The  dynamic  effects  of  orbital 
eccentricity  and  the  earth's  oblateness  were  neglected. 

In  the  dynamic  simulation  using  Eq.  (.87)  Bainum  and 
Kumar  used  a  commanded  length,  L,:,  exponentially  increasing 
with  time  as  m  Kalaghan  et  al .  [13]  .  with  suitable  modifi¬ 
cation  for  the  deployment  ,  i.e. 

Lc  =  10  *  S4S?0  (  l-e-t/-° )  (88) 

where  p  is  a  positive  constant.  With  this  commanded  length 
used  in  conjunction  with  the  optima-  feedfcacJ;  gains,  deploy¬ 
ment  requires  a  moderate  duration  :f  time  and  has  better 
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eristics,  w.iile  comparing  tnis  result  with 
[Baber  et  al . ] 1 and  BAG  'Kalaghan  et  al . ] 
dep -oyment  simulation  without  aerodynamics 
- ane  mot  i  in  required  appr animate ly  G . n  hr.  to 
ateiiite  to  IJOhm,  wnereas  the  SAG  deployment 
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using  Rupp's  [10]  law  with  atmospheric  effects  included 

required  approximately  10  hr.,  to  deploy  the  subsatellite  to  100  km. 

3.2  Stationkeeping 

This  task  involves  the  maintenance  of  the  subsatellite 
in  the  desired  equilibrium  configuration.  For  the  present 
proposed  aerodynamic  test  mission,  the  desired  equilibrium 
configuration  involves  a  sutsateilite  deployed  IDO  km  below 
(or  above;  the  Snuttie  Grbiter  (Shuttle  Orbitsr  altitude  - 
220  km;  along  the  local  vertical.  Under  the  influence  of 
external  disturbances,  the  subsatellite  deviates  from  this 
desired  position. 

Fainum  and  Kumar  [15]  used  optimal  control  law  theory  to 
investigate  station  Keeping  for  the  Shuttle  Tethered  Sub- 
satellite  system.  The  system  response  to  various  initial 
'cnditicns  with  all  the  external  disturbing  forces  absent 
was  studied.  Typical  results  shewed  that  the  time  constant 
for  tile  least  damped  mode  corresponding  to  the  optimal  gains 
0.93  •  rbit )  was  shorter  than  "hose  associated  with  Rupp's 
i aw  1.351  orbits,  in-plane  tuning:.  In  Rupp ' s  control  law, 
the  length  and  length  rate  gains  were  selected  for  m-plane 
-r.d  ■  'it- :  f-pl  ate  control  based  on  the  system  natural  fre- 
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U  =  -  ( R_1BtP)X  =  -  KX  {  90) 

The  system  response  with  lower  overshoots,  shorter  settling 
limes,  and  with  comparable  power  and  tension  levels  was 
obtained  by  applying  modern  control  theory.  Maximum  tensile 
acceleration  for  this  case  was  0.525  m/sec-  for  the  optimal 
control  and  0.505  m/sec-  for  Rupp's  [10]  in-plane  tuning 
control . 

When  the  effect  of  atmosphere  on  the  tether  was 
considered  with  Rupp's  [10]  control,  the  system  response  had 
Less  damping  than  with  the  optimal  control.  This  reflected 
the  greater  stiffness  in  the  optimally  controlled  system. 
Also,  Rupp's  control  had  a  tendency  to  pull  the  subsatellite 
to  a  higher  altitude,  unlike  the  optimal  control  which  had  a 
tendency  to  deploy  the  subsatellite  further  into  the  atmc- 


s ohe re . 


is  evident  that  the  control  law  based  on  the 


Linear  regulator  theory  resulted  in  a  superior  performance 
when  compared  to  Rupp ’ s  control  law  m  the  station  keeping 


.HOC?  . 


“hen  the  previous  work  was  later  extended  to  Platform 
’ether  ed  -  Sub  satellite  systems  by  Eamum,  Woodard,  and 


Jvang  [  *  ] ,  where  a  "-dimension  mathema" xa. 


.ei  -t  open 


and  closed  loop  in -orb it  plane  dynamics  of  a  space  Pi.atform- 
?e  iherec-Subsatellite  system  was  treated.  The  control  was 
-earned  tc  be  provided  by  both  modulation  of  the  tether 


n&v:*:.-  vn:;  ca.atf  :rm 


pitch  angle,  respectively,  were, 

Tr  =  K  e  +  K  ,e '  +  K  a  +  K  ,a'  +  K  .0  +  K  Q,0' 

u  ee  ee  ea  v  ea1  v  e9  v  e0  v 


91 


=  Ka  £  +  K.  ,£’  +  Kq  a  +  Ka  ,a 1  +  K  0  +  K  ,0’  (  92 

9  0£  ftp1  firv  tt  fin1  T7  AA  Tr  Qfi'  tt  x 


0£ 


9a  v  0a'  v  60  v  90 '  v 


where  £  =  L/L,~  -  l,e'  ,  a  =  a-a  ,  a',  9  =0-0  .0' 

v  eq  v'  v  eq'  v 

and  K 1 s  are  state  variables  and  gains,  respectively.  The 
angle,  9  ,  describes  the  orientation  of  the  platform  with 
respect  to  the  local  vertical  and  the  angle,  a,  represents 
the  angular  displacements  of  the  tether  line  relative  to  a 
local  normal  in  the  platform,  Lc  is  the  nominal  reference 
length,  and  the  subscript  "eg”  refers  to  equilibrium 
values . 

The  numerical  results  showed  that  tether  line  swing 
motion  was  damped,  requiring  about  1.75  hr  to  reach  the 
nominal  value,  whereas  the  platform  pitch  motion  was  damped 
out  within  approximately  1.0  hr  when  corresponding  feedback 
gains  and  initial  conditions  were  determined. 

It  was  proved  that:  (1)  within  the  linear  range  the 
system  is  controllable  with  momentum- type  control  on  the 
platform  and  with  tension  modulation  m  the  tether  line;  ( 2 ■ 
the  linear  system  is  observable  with  tether  length  and 
.Leagtu  rate  measurements  only;  i  Z)  the  tether  attachment 
point  offset  increases  the  system’s  natural,  coupling  and 
improves  transient  performance  m  the  least  damped  mode,  but 
at  the  cost  of  slightly  larger  control  force  amplitudes. 

Continuing .  Far.  Ruyiag  and  Sair.um  |5  ]  developed  a  3- 


dimensional  mathematical  model  of  the  open  and  : hosed  loot 


3  c  ace  Tethered.-?  1st  form-: 


:em 
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and  used  the  same  fcrm  of  the  control  law  as  in  Ref.  [15]. 

It  was  assumed  that  the  control  could  be  realized 
through  appropriate  modulation  of  the  tension  in  the  tether 
line  and  the  momentum  type  controller  for  the  platform 
pitch,  roll  and  yaw  rotation,,  i.e. 
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( m -plane ) 
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where  the  angles  0  ,  <f>,  ip  ,  and  0',  <f>',  \p'  are  platform  pitch, 
roll  and  yaw  angles  and  their  rates,  respectively,  <4,a'  ,  y  ,y' 
are  tether  line  m-plane  and  out-of-plane  swing  angles  and 
angular  rates,  respectively. 

The  numerical  results  showed  that  the  platform  pitch 
angle  and  tether  line  swing  angle  damped,  both  requiring 
about  ..5  hr  for  the  initial  conditions  selected.  For  out- 
of-plane  motion  the  platform  roll  and  yaw  angles,  and  the 
tether  line  swing  angle  damped  requiring  about  1.5  hr,  3  hr 


and  _.5  hr,  respectively. 

For  the  case  where  the  tether  attachment  point  offset 
wa.  only  a .eng  the  roll  an is ,  it  was  verified  that  both  the 
m-plane  and  out-of-plane  sub satellite  systems  are 
:: r.t roll  able  wren  tension  modulation  or.  the  tether  and 
-.unent  un  ty?*-  control  are  available.  Both  subsystems  are 
-r  su_-  _f  the  length  of  “he  tether,  and  the  platform 
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rotation  angles  together  wich  their  rates  are  available. 

The  terher  attachment  point  offset,  which  is  the  source 
of  the  system's  natural  coupling,  is  an  important  factor  in 
establishing  system  controllability  and  observability.  For 
the  case  of  no  attachment  offset,  the  rotation  of  the 
platform  will  not  affect  the  subsatellite  cut-of-plane 
swing;  in  other  words,  the  effect  of  higher  order  terms 
should  be  considered,  or  other  means  of  control,  such  as  by 
placing  an  actuator  on  the  subsatellite  to  control  the 
tether  line  cut-cf-plane  swing,  should  be  augmented. 

The  investigation  of  the  effect  of  tether  flexibility 
on  the  : in-plane)  stability  regions  as  a  function  of  the 
■ether  tension  control  parameters  during  the  station  keeping 
was  further  developed  by  Liu  Liangdong  and  Bainum  [8  ]  where 
an  alternate  opcimal  control  strategy  which  included 
additional  feedback  of  the  first  vibrational  mode  and  its 
rate  was  introduced.  The  formulation  of  tension  level 
control  was  m  the  form 

T  =  -  (K  £  +  K  ,e'  +  K  a  +  K  ,a'  +  K  n,  +  K  A1)  (96  ) 

e  e  a  a'  1  1 

where  is  the  first  flexible  modal  amplitude  (non- 

dimeus icaal ;  stare  variable  and  its  rate,  respectively.  The 
"'/stem  is  control _abi-  using  only  tether  tension  control  and 
-.Iso  Dbser"able  with  the  measurement  .  f  s,e •  a, a’,  or  only 
£ / £ 1 ;  n]_'n[  tre  available  either  through  estimation  or  by 

dir act  me  - surerent . 

A  .-cording  to  n  imeri  :al  comparison  of  che  transient 
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responses  it  was  seen  that  Rupp 1 s  control  law  could  be  used 
to  control  the  the  in-plane  swing  angle  succec siui 1 y  during 
stat.cn  keeping  but  it  was  not  very  effective  for  damping 
the  in-?lane  vibrations.  The  transient  responses  for  Bainum 
and  Kumar's  !  15]  optimal  law  based  on  t,e'  ct,  a'  were 
faster  than  those  for  Rupp’s  10]  control  law.  The  further 
improvements  m  transient  response  in  both  the  m-plane 
swing  angle  and  the  first  vibration  modal  amplitude  were 
made  by  including  the  state  feedback  of  the  first 
vibrational  mc-de  into  the  optimal  control  law.  An  impove- 
ment  was  also  apparent  in  the  damping  of  the  second  mode  due 
to  the  coupling  between  "lie  first  and  second  modes. 


3.3  Retrieval 

The  retrieval  is  basically  an  unstable  procedure. 
Retrieval  can  be  carried  cut  by  letting  the  commanded 
length,  Lc ,  reduce  with  time.  bc  can  ce  decreased  in  steps 
[Rr.ppj-  or  it  can  be  an  exponentially  decreasing  continuous 
function  of  time,  such  as 


Li,-  =  L.-.e-t/P  .97) 

In  either  case,  the  rotational  as  well  as  vibrational 


in-  are  un. 


m  the  absence  of  active  control, 


-cause  to-  negative  damping  inf  reduced 
\  L/L  and,  m  practice ,  tl 


during  retrieval  is 
e  camping  level 
rev" ired  t  c  guarantee  stability  is  not  always  available. 
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aw-;  were  used  i.s  retrieval  but.  large  amplitude 
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45  respectively,  resulted,  depending  on  the  initial 

With  3  si  nun  and  Kumar's  [15]  optimal  control  strategy 
for  retrieval,  the  in -plane  re-'t  :nse  is  beater  than  that 
obtained  by  Baker  et  al.  fl2] ,  however,  cr  -  ;f-r  ae  motion 
i s  hardly  improved . 

To  restrict  out-of-plane  oscillations  to  reasonable 
bounds,  nonlinear  control  strategies  must  be  used.  Xu  et 
al .  [  16]  showed  that  a  satisfactory  length  rate  law  is  of 
tre  form, 


=  K  C 1  +  K  ,c: '  +  K  ,  y ' 2  ] 

T .  Q  '  y  ' 


■98 


L  *  9 

wrers  K^,  and  K  ,  are  negative  constants  K0  is  a  negative 
function  of  the  true  anomaly  9.  a  is  the  n.-plane  tether 
swing  angle,  and  y  the  out-of-plane  tether  swing  angle. 
Modi  et  al .  [17]  considered  some  nonlinear  control 

strategies  with  a  tension  control  law  of  the  form, 
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T  =  K^L  +  F^.L'  +  K  ,y'2  +  T 
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and  desired  final  condition  together  with  a  cost  function 
that  penalizes  in-plane  and  out-of-plane  deviations  of  the 
motion  and  uses  the  commanded  tether  length  as  the  means  of 
control . 

It  is  sought  to  minimize  a  cost  function  i.e. 

fcf 

J=  b[x(tf),t  ]  +  /  {  L[x(t)  ,u  (t)  ,t]}dt 

o  (100  ) 

'ey  appropriate  choice  of  the  scalar  control,  u(t),  where  t0 
and  if  are  given.  4>  [>:(tf),t«]  is  the  terminal  cost  and 
L[x;t),  u(t-,tj  is  the  Lagrangian.  These  are  most  common  of 
the  quadratic  forms: 

[  x ( t  f  )  ,  tf]  =S  Wf  {  xf-x^)  j_-  (  101 ) 

1  [  x  (  t )  ,  U  (  t }  ,  t  ]  =  I  c  1xi2  (  t )  +  F;  U2  (  t )  (  ) 

where  Wj_ ,  Qx  and  R  are  weighting  factors  that  indicate  the 
importance  of  minimizing  the  associated  state  component,  x-j. , 
or  control,  u.  and  xa  is  a  desired  final  state  such  as  the 
final  desired  tetl.er  length.  With  the  implementation  of  a 
first  order  conjugate  gradient  method  the  optimal  L  control 
was  obtained.  The  numerical  simulation  results  showed  that 
by  suitable  choice  of  state  and  control  weightings  the 
i  irrational  motion  could  be  significantly  reduced,  and  in- 
pi  ane  ; -  ciilat:  ns  ■••ere  mere  readily  attenuated  than  those 
f:r  the  out-of-plene  motion. 

Tie  more  realistic  model  was  net  considered  in  this 
irves _ i gat i  :r.  suer  as  including  one  "ransverse  motion  of  the 
t-tr.e.  ar. :  t..e  e freer  .  f  atmispnerir  forces  arc  eccentricity 


45 


cf  the  orbit. 


In  order  to  further  improve  the  performance  of  the 
system  an  additional  nonlinear  tension  control  law  was 
introduced  by  Liu  Liangdcng  and  Bainum  [  8]  which  is  of  the 
f  o  rm : 

T  -  -  (  m3r;H|-  )a)  -  [  F  £-*-FKj  AL-FI^hL'  +  L(FK3Ct~ 

+  FK4Y2  +  FK5  y'  2  )  ]  (  103  ) 

where:  sn3 ,  mt  are  the  mass  of  subsatellite  and  tether, 

respectively;  ^  is  the  orbital  frequency;  AL  and  iL1  are 
the  difference  between  the  tether  length  and  some  reference 
length  (and  its  rate;;  the  angles  ot,Y,a,Y'  are  tether  in- 
plane  and  cut-cf-plane  swing  angles  and  angular  rates, 
respectively;  and  Fu^  are  optimal  control  law  gains. 

Because  it  was  difficult  to  use  strictly  analytical 
methods  to  der.ve  centre i  gams  for  such  nonlinear 
equations,  the  cost  function  was  selected  as  oelcw: 


fQ  iC-  +  Y^)dt  ;/Tf 


(  104 ) 
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compared  with  Rupp's  [10]  control  law,  and  the  transient 
responses  were  also  better  than  those  based  on  the  control 
law  in  Ref .  [ 17 ] . 

The  tension  control  is  unreliable  during  the  terminal 
stage  of  retrieval  when  the  tension  becomes  very  small 
because  of  the  small  length  of  the  tether.  The  tension 
might  even  become  zero  (slack  tether)  due  to  the  longitu¬ 
dinal  oscillations .  To  overcome  this  difficulty  Eanerjee 
and  Kane  [19  ]  proposed  that  natural  tether  tension  could  be 
augmented  with  satellite-based,  tether-aligned  thrusters  and 
tr.at  these  thrusters  would  be  capable  of  stabiizing  and 
speeding  up  tne  retrieval  process.  The  thruster  augmented 
torque  control  was  of  the  form. 

T.- -T~  =  K  a  +  K  a'  +  X  y  +  K  y’  (107) 

^  '-a  a  a  y  Y 

where  Tc  is  the  torque  proportional  to  a' and  y'  . 

A  summary  of  the  control  laws  used  by  the  various 
lii v ^ s t i crs 1 1 oil o  is  giv^n  in  3 

3.4  Recommendation  Remarks 

(1)  Because  mcst  cf  the  useful  missions  are  cairied  cut 
during  the  station  keeping  phase  for  Shuttle/Platform 
Tethered  Bubsatei Lite  systems,  the  investigation  of  control 
laws  shot. _d  be  foci 
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augmentation  of  the  tether  tension  control),  it  may  still  be 
difficult  to  implement  an  efficient  and  successful  retrieval 
for  the  tethered  reflector  system. 

( 2 )  For  deployment,  the  active  tension  modulation  schemes 
proposed  by  Rupp  ft-0  J  and  subsequent  investigators  are  more 
efficient  than  the  purely  passive  scheme  advocated  by  Kane 
and  Levinson. 

(2)  Out  of  the  various  active  control  laws  for  proposed 
orbiting  tethered  reflector  systems,  Rupp's  [10]  control  lav/ 
is  the  most  basic  and  is  effective  in  controlling  in-plane 
motion,  but  not  adequate  for  out-of-plane  motion  control. 

( 4*.  Kissel's  ’Baker  et  a.1 .  ]  tl2]  tension  control  law,  based 
on  a  :ombination  of  exponential  and  uniform  rate  of  tether 
length,  as  a  commanded  length  rate  can  be  used  both  for  in¬ 
plane  and  out-of-plane  motion  control  during  deployment.  It 
was  seen  that  the  duration  of  time  for  m-piane  deployment 
was  reduced  as  compared  with  that  in  Rupp  [10]  ,  and  damping 
characteristics  for  both  the  case  of  station  keeping  and 
deployment  were  better  than  those  for  Rupp. 

(5)  Eamum  and  Kumar's  [15]  control  lav;  based  upon  the 
linear  regulator  problem  of  optimal  control  theory  is 
suitable  f:r  adaptation  with  proposed  Shuttle  -Tethered 
systems,  ..tore  specif icaiiy : 

For  station  :eeping  purposes  at  altitudes  where 
atmospheric  effects  ere  negligible,  control  laws  result  in 
si  ant  response  to  initial  perturbations  as 
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Kissel’s].  The  tether  tension  and  power  levels  required  for 
such  control  do  net  exceed  previous  requirements  (for  the 
TSS  system':  . 

For  steady  state  station  keeping  requirements  where 
atmospheric  effects  and  eccentricity  may  be  important,  the 
optimal  control  law  with  gains  based  on  optimal  control 
theory  can  be  used  to  bring  the  system  to  an  in-plane 
equilibrium  position  and  tension  level  which  reflects  a 
balance  between  the  gravity-gradient  and  aerodynamics 
torques  (forces). 

For  moderate  duration  deployment  the  same  form  of  the 
optimal  law,  where  the  actual  gains  are  adapted  continuously 
to  the  commanded  length,  can  result  in  improved  damping 
(settling)  characteristics  with  small  amplitude  initial 
excursions  m  the  m-plane  swing  angle. 

if;  The  same  form  of  the  optimal  control  law  as  that  in  [15] 
can  be  used  extensively  in  both  2-dimer.sional  dynamic  models 
[Bainum  and  Woodard]  14  ]  and  3-dimensionai  dynamic  models 
(Fan  Ruying  and  Bainumj!i5],  of  the  Platf orm-Tethered-Sub- 
sateliite  systems  during  stationkeeping. 

{“}  when  considering  tetr.er  mass,  tether  flexibility,  aero- 
ryramic  force  on  the  tether,  and  eccentricity  of  the  orbit, 
i.e.  a  mere  complex  dynamic  model  for  tethered  systems,  Li.i 
liangdcng  and  Bainum'. a  alternate  optimal  control  lav;  made 


farther  improvements  in  the  trails 


respor.se 
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To  sum  up,  control  laws  based  upon  optimal  control 
theory  offer  the  greatest  potential  for  applications 
involving  proposed  orbiting  tethered  reflector  systems. 
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Summary  of  Tether  Control  Laws 
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Investigators  Motion  Type  Form  of  Control  Law 

controlled  Control  _ _ 

Bainum  2d.  rota-  optimized  TL-K££6  +K£e*6/+I<^v+^ 

Woodard  tions  and  tension, 

Juang  stretch  momentum  +K6£'£/+Ket/^^Keo(;^Kd^v0+Kd^/d 
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Investigators  Motion  Type  Form  of  Control  Law 
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4.  NONLINEAR  DYNAMIC  EQUATIONS  OE  A  IE. THERE!) 

ANTENNA/REFLECTOR  IN  ORBIT 

4.1  Introduction 

The  linearized  equations  of  the  orbiting  tethered  antenna  system  have  been  ob¬ 
tained  in  Chap.  II,  and  active  tether  control  laws  (I ,QR)  based  on  the  linear  model  have 
been  efficiently  used  to  apply  active  control  during  stationkeeping  of  the  system. 
However,  for  the  deployment  abd  retrieval,  the  linear  system  model  may  not  represent 
the  physical  situation  accurately  any  more  and  the  active  control  laws  based  on  this 
model  certainly  may  not  be  as  effective  during  deployment  and  retrieval  as  during 
stationkeeping.  This  is  due  to  the  large  slewing  angles  and  inherent  instability  of  the  out- 
of-plane  motion  of  the  tether.  Furthermore,  second  order  terms  in  length  rate  are 
directly  coupled  with  out-of-plane  modal  amplitude  terms.  Hence  to  damp  the  out-of¬ 
plane  motion  using  length  rate  control  (  tension  control  )  and  to  simulate  the  dynamic 
behavior  of  the  system  during  the  deployment  and  retrieval,  it  is  necessary  to  use  the 
nonlinear  equations. 

The  general  dynamics  of  a  tethered  system  is  rnthei  complex  and  hence,  early 
dynamical  models  were  based  on  a  number  of  simpllving  assumptions.  An  overview  of 
the  development  in  this  area,  parlicuraly  system  models  and  proposed  control  laws,  has 
been  given  by  Misra  and  Modi  (14]  and  Bainum  am1  Kunw  |  15],  'Hie  syslcm  models 
have  grown  from  initial  massless  idealized  tether  models  m  vomplex  representations  en¬ 
compassing  all  the  tether  vibrations  (  flexibility  )  and  cm!  hnd\  motions,  as  exemplied  in 
the  model  by  Misra  and  Modi  f  1 4 ] . 

As  for  the  iethered  antenna/rcflector  system,  die  translational  motions  of  the 
subsatellite  and  transverse  vibrations  of  the  tether  will  affert  'be  rigid  body  motion  of  the 
orbiting  antenna;  therefore,  all  these  effects  will  be  included  m  'he  formulation  of  the 
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nonlinear  system  equations  for  the  further  study  of  simulating  t he  system  dynamic  be 


havior  and  applying  active  control  during  the  deployment  and  retrieval. 

4.2  The  Assumptions 

The  following  assumptions  are  made  to  develop  the  model  equations: 

a) .  The  antenna  mass  is  far  greater  than  that  of  the  tether  nnd  the  subsatellite; 

consequently,  the  center  of  mass  of  the  system  mnv  he  taken  to  coincide  with 
that  of  the  antenna. 

b) .  The  tether  is  assumed  uniform  with  a  constant  mass  per  length. 

c) .  Longitudinal  stretching  is  not  considered  .  and  longitudinal  vibration  is 

neglected  compared  with  the  transverse  vibrations. 

d) .  The  shell  is  considered  to  be  a  rigid  body. 

e) .  The  subsatellite  is  considered  to  be  a  point  mass. 

f) .  No  random  inputs  or  unknown  disturbances  are  considered. 

g) .  Only  first  order  gravity-gradient  effects  are  considered  and  the  orbit  is 

assumed  circular. 

4.3  Kinematics  of  the  System 

The  coordinate  systems  used  in  the  development  of  the  system  equations  of  motion 

are  shown  in  Fig.  7.  Op  Xo  Yo  Zo  is  an  orbit-reference  centered  at  the  center  of  mass  of 

the  shell,  Op  ,  with  Op  Xo  along  the  local  vertical  and  Or  >  along  the  orbital  angular 

velocity  direction,  Op  Zo  along  the  orbital  tangent  velocity  hhcc'ion. 

OpXpYpZp  is  a  shell  body  reference  frame,  RP  .  ''  here  Or  XP  ,  ()PYP  ,  Op  ZP  are 

principal  axes  of  the  shell.  OXt  Yt  Zi  is  the  subsatell  te  undeformed  tether  reference 

frame.  Rt  ,  with  OXt  along  the  undeformed  tether  line,  "  lice  O  is  the  point  from  which 

the  tether  is  deploying  or  retrieving.  The  coordinates  of  O  i-i  the  shell  frame  are  (h,  0, 

0). 

The  Euler  angles  '1/,  0,  rt>  are  the  yaw.  pitch  and  roll  males  of  the  shell,  respec- 
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tively.  a,  y  are  pitch  (in-plane)  and  yaw  (  out-of-planc )  angles  of  the  tether. 

For  convenience,  the  transverse  vibrations  of  the  tether  nre  expanded  in  terms  of  a 
set  of  admissible  functions. 

V=£$n(x)Bn(t)  W=£4>n(x)Bn(l)  (108) 

where  4>n(x)=sin(  ntix/L)  ,  v — out-of-planc  displacement  of  the  tether,  w — in-plane 
displacement. 

Therefore,  the  whole  system  tias  the  following  degrees  of  freedom: 

vi/,  0,  cb  - rigid  body  motion  of  the  shell. 

a,  y  - translational  motion  of  the  subsateHhe. 

L  - length  of  the  tether. 

Bn.  Cn  - transverse  vibrations  of  the  flexible  'ether. 

The  transformation  matrices  from  Or  Xr  YP  Zr  to  n;.  X  •  Y..  Zo  and  C)  XtYiZt  to 
Op  XP  YP  ZP  are  given  by 


O)  Gb  -  CO  S6 

MP  =  -S4/S0Gt>+ Ql/S4>  Srl/SOSb  +  CVrOb 

QlrSe  QM-  S^l/Scb  -  QlrSO  S<M-  StlrG't 


sn 

Stir  CO  (109) 

ntrco 


M 


C  a  C  y  -CaSy  -Sn 

S  y  Cy  0 

S  a  C  y  -  SaSy  C  nr 


(110) 


(Hi) 
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X 

“  h  ‘ 

X 

p 

t 

y„ 

— 

o 

+  M 

t 

y, 

z 

L  r  J 

0 

z 

_  t  _ 

(112) 


where  C  —  cos,  S  --  sin 

The  angular  velocity  of  Ihe  shell  is  given  by 

/y  /y  S' 

cop  =  12  x  ip  +  12  y  jp  +  12  7  kp  (113) 

where 

12  x  =  <1/  C  0  C  4>  +  0  S  <!>  +  (i>c  (  -  n  iif  S  a  C'  o  t  <  i’i  S  (li  } 

12  y  =  6  C  <b  -  >Ir  C  0  S  <h  f  wc  t  S  U«  SOS  S  ( ’  >'■  C  <b  I  (114) 

12  *  =  < t>  b  Ur  SO-  coc  S  \!r  C  0 
The  angular  velocity  of  the  tenter  is  given  by 

A  ,A  A 

(Ot  =  (Ox  it  *t~  OJy  jt  +  (Ox  kt  (115) 

(Ox  =  Ca  Cy  12  x  +  Sy  12  y  f  Sa  Cy  12  •  i  ’ 
coy  =  -  Ca  Sy  12  x  +  Cy  12  y  -  Soc  Sy  12  /  t  <''/  ,v  (1 16) 

coz  =  -  Sa  12  x  +  Ca  12  z  +  y 
4.4  Dynamic  equations  of  the  System 
4.4.1  Rigid  Body  Motion  of  the  Shell 

The  Euler-Newtonian  method  is  used  lo  develop  he  d\n  *»iii<:nl  eouations  of  Ihe  shell 
motion. 

The  angular  momentum  of  the  shell  is  given  hv 

N  =  1  xx  I2x  ip  t  lyy  12yjp  s  \77il7  kp  (11  ) 

The  time  derivative  of  the  N  is  written  as 

N  =  f  Ixx  1 2x  "  (  (vv  'I77  )  12y  1  2z  |  ip  u  (  lyy  12y“ t  I  '7-  I  <x  1 1 1  x  1  2  ■  I  I  !  1  .•  ’  1 2.’  "  (  I  xx  ~lyy )  1 2x  1 2y  )  kp 

There  are  two  forces  acting  on  the  shell,  one  '■  the  in  'national  force,  the  other 


(nsj 


tether  tension  force.  T.  Therefore,  the  torque  exerted  <>n  th"  '•'hell  i« 


(119) 


L  =  Lg  +  Lt 

where  Lg  is  the  gravitational  torque  and  Lr  is  the  'orque  ot  die  tether  tension 
force. 


L r  =  -hT sin«  cosy  ]r  +  hTsin-y  k,.  ( 1 20) 

[(ere,  vve  only  considci  first-order  gravitational  force*--  It  is  well  known  that  the 
first-order  gravitational  force  acting  on  any  point  Pfx.y./l  in  'lie  sficll  can  be  expressed 
[7  ]  in  the  shell  frame  as: 


’  3  CD  Gb  -1 

-  3  CO  S6  Gb 

-  1  so  fTt  f  V, 

A 

X 

elf  =  ox-dm 

-3  CR2S6  Gb 

3  CDS6-  1 

t  SO  ( Tt  S,t, 

/\ 

y 

-3  SO  CD  C<b 

3  SO  CO  Sfb 

a  sn  t 

a 

7 

(121) 


Hence,  the  torque  of  the  gravitational  force 


lr,=f  rxdf  = 


3  SO  on  S6  (  1  i  1 

yy 

-3  SO  CO  Gt>  (  I  -  I  i 

^  77  xx 

-3  CO  Srt  Ob  (  I  -  I  i 


(1221 


According  to  the  Euler-Newtonian  equation,  we  obtain 


i  r.  -  (  i  -  i  )  a  a  =  3  <o 2  so  cti  s.' 

xx  x  yy  7.7  y  2 


il  -  (  I  -  I  )  a  a  =  -  ItTSa  (Jv  t  .>sn  Ml  o'.  (  I  -  I  ) 

y  77  XX  X  2  '  77  XX  (  kJ) 


1  U  -  (  I  -  I  )  a  a  =  Il  l'S  y  ~  3  ..)  nt  <v,  ■ 

77  1  XX  VY  X  V 


The  above  arc  the  nonlinear  dynamical  equations  |m  'he  -!■  >l| 


4.4.2  Translational  Motion  of  the  Subsystem  (tclltc  uul  s'l^mellhe) 

The  La n grange  approach  is  used  to  develop  the  dv naurmti  equations  of  the  transla 
t  tonal  motions  of  the  subset  cm.  Thercfo—  we  need  ■<»  cacui  ue  he  kinetic  and  potential 
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energy  of  the  tether  and  the  subsatellite. 

A).  Kinetic  Energy 

The  velocity  of  any  point  Q(x,y,z)  on  the  tether  is  of  the  following  form: 

__  _  • 

v  =  Vo  +  w«  x  n  +  7.  t  124) 

where 


V 


-So.  Cy  H  +  Sy  Cl 

i. 

_  a 

i  _ 

y  7. 

A 

1 

i 

So.  Sy  il  +  Cy  Cl 
y  7. 

+ 

0 

} 

J 

A 

-  CaH 

L  y 

.  0  . 

k 

_  '  - 

(125) 


w 


x  r  =  < 

'  1  l 


A  _ 

0 

w 

-v 

(0 

i 

:< 

A 

-w 

0 

X 

(0 

i  t 

A 

V 

~x 

0 

/ 

l* 

_ 

0) 

(  126) 


r  = 

i 


T 

0 

i 

t 

A 

V 

j, 

A 

k 

w 

t 

t  127) 


Hence,  the  kinetic  energy  of  the  tether  is  given  b\ 

T  ,  =  ]  p  M2^x  =  /  p  flfdx  +  4  f  p  ( to  xr  |  ih  '  4  /  p  i 7  1 2dx 

n  o  "*  . . 

i.  _  .  i.  ,  i.  _  _  .  (128) 

+  /  p  Vo- (to  xr  )dx  +  /  p  Vo-r  dx  +/  p  (w  xr  i-r  dx 


where 


T>  =4/PN‘fdx 


m  2  2  2  2  .  j. 

oJ  f  h  n  +h!l  t!.+  2h  I  '  mrv  i> 


-  <h-.r  cosy  12  )  ]  ( 129) 
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T  =  -|-/p  pit  xFt  I  dx  =  OJ  (r  ; p 


v2+v^  -X'  -xiv 

-xv  x2t  »•  -'-w 

-XVV  -V"  V2 


dx  co 


m  1  2  2  ,  I2  I 

1  r/"  n  «  r~-  \  .  A  «  (  i  1 


L2  ,  1  T-  n  \.  .2 


=  ~T  1  4  1  T-  +  Ts  n„K 


+  2_L  Z_L  D  B  w  w  +  2  c  co  co  -SBC  -o  co 

TT  11  nxy  TT  n  n  x  /  nv 


(130) 


T  =  -j  I  p(  |v  |2+  |w|2)dx 


where  v  =1  sirvT-B  -  I  n  x  cox  B 

L  n  j2  1 


(131) 


~  rmx  TT  L  ^  M  „  ITT7X  f  • 
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T  s  =  /  pVo  .  r (  dx  =  — Hy  [  sin  a  siny  S— R(1  c<m  z'-tLjjT  Cn  + 
sinasirry  (  J3n  -  cosot  (  -}■-  )  S-'-ti"  1  f>+ '  C  n  ]  (133) 

+  C0SY  nz  f  nn  +  (  n  | 


T  =  J  p  (  £3  x  r  )  .  r  dx  =  to  T  f  p 

6  J  v  t  t  t  t  * 
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(134) 


The  velocity  of  the  subsatellite 


V  =  Jhll 
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(135) 


Hence,  the  kinetic  energy  of  the  subsatellite 


2  2  22  .2 


T  =  — i-  m  |V  |2= — i —  m  [h  fl  +  h  11  f  I.  +■  I  to  *  |  2h  «in nt  cosy  11  I, 

<  Z  <  1  <r  1  Z  *:  1  y  7.  '  •  y 


(136) 


+  2hsinyll  L  +  2hLll  (cosot  11  +sintt<in"v  <•>  i*  ’hi  <-os-v  ll  oj 

t.  y  y  '  ' 


13),  Plastic  Potential  Hnergy  of  the  Tether 

Suppose  the  tether  is  an  Fule:-I3crb<nri  beam  I  bore  fore,  the  elastic  potential 
energy  is  of  the  fol'o’’  ing  form: 
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E  =  -4-J  El  (  V2  +  w2  )  dx  = 
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C).  First-order  Gravity  Gradient  Field 

Since  the  origin  of  the  orbital  reference  frame  move*  along  a  free-fall  trajectory, 
the  only  gravitational  forces  acting  on  the  subsystem  arise  from  the  gravity-gradient  field. 
The  gravity-gradient  force  terms  are  obtained  by  a  Taylor- series  expansion  of  the  gravity 
field  about  the  free-fall  trajectory.  The  first-order  terms  ol  ibis  series  are  well  known. 
Since  Bn  /  L  <tc  1  ,  Cn  /  L  1  ,  the  effect  of  the  transverse  Gbration  of  the  tether  on 
the  gravity  force  is  neglected.  Applying  these  terms  to  a  mass  particle  of  si?c  dm  results  in 
the  following: 
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c  J 

dF 

=  -  a>2  z 
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Z  C 


where  (  x,  y,  7.  )  are  the  coordinates  of  the  particle  dm  in  the  orbital  reference  frame,  toc 
is  the  orbital  angular  velocity. 

Summing  up  the  forces  over  all  mass  particles  of  the  dvnamr.nl  system  yields  t he  first— or 
der  gravity-gradient  terms  as: 
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Since  the  antenna  mass  is  far  greater  than  that  of  the  tether  and  subsatellitc,  it  is 
assumed  that  the  Euler  angles  of  rigid  body  motions  of  'lie  satellite  and  their  rates,  the 
vibrating  mode  shape  and  their  rates  are  small  (  first  order  terms  ).  that  is  ,  rb,  0,  <b,  Q  *  .  Ay, A 

,  C  ,  .  fl.i  ,  Cn  <  1  .  Omitting  the  third  order  terms  and  above,  we  obtain  the 
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following  dynamical  equations: 
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4.4.3  Mode  liquations  of  the  Tether 

The  vibrating  modes  are  coupled  with  (lie  other  degrees  of  freedom.  They  satisfy 
the  following  equations: 

Out-of-plane  mode 
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In-plane  mode 
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4.5  Conclusion 

The  nonlinear  dynamical  equations  of  the  tethered  amentia  system  have  been 
obtained.  We  can  see  that  all  the  degrees  of  freedom  are  coupled  in  the  equations.  From 
Eq.  (142)  and  Eq.  (143),  it  is  seen  that  the  in-plane  and  mit-of-plane  motions  are 
coupled  through  second-order,  and  also  coupled  with  the  flexibility  of  the  tether.  The 
dynamical  behavior  of  such  a  complex  system  (  including  altitude  motions  of  the  satellite 
and  flexibility  of  the  tether  )  has  never  been  studied  before.  Next  step  of  our  research  will 
concentrate  on  the  simulation  of  the  dynamical  heinr  iot  of  'he  system  during  the  deploy 


ment  and  retrieval. 


Chapter  4 


Some  integrals  used  in  this  development 
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5.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  system  linear  equations  for  the  motion  of  a  tethered 
shallow  spherical  shell  in  orbit  with  its  symmetry  axis 
nominally  following  the  local  vertical  are  developed.  The 
shell  roll,  yaw  tether  out-of-plane  swing  motion  and  elastic 
vibrations  are  decoupled  from  the  shell  and  tether  in-plane 
pitch  motions  and  elastic  vibrations.  The  neutral  gravity 
stability  conditions  for  the  special  case  of  a  constant 
length  rigid  tether  are  given  for  in-plane  motion  and  out- 
of-plane  motion.  It  is  proved  that  the  in-plane  motion  of 
the  system  could  be  asymptotically  stable  based  on  Rupp's 
tension  control  law,  for  a  variable  length  tether.  However, 
the  sytem  simulation  results  indicate  that  the  transient 
responses  can  be  improved  significantly,  especially  for  the 
damping  of  the  tether  and  shell  pitch  motion,  by  an  optimal 
feedback  control  law  for  the  rigid  variable  length  tether 
model.  It  is  also  seen  that  the  system  could  be  unstable 
when  the  effect  of  tether  flexibility  is  included  if  the 
control  gains  are  not  chosen  carefully.  The  transient 
responses  for  three  different  tension  control  laws  are 
compared  during  typical  station  keeping  operations.  The 
transient  responses  can  be  further  improved  by  including  the 
state  feedback  of  the  tether  vibrational  modes  into  the 
optimal  control  law,  especially  for  the  damping  of  the 
tether  vibrations. 
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Extensions  to  the  present  study  could  consider  the  effect  of 
the  shell  flexibility  on  the  system  stability  and  control 
and  some  kind  of  active  control  could  be  introduced  (in 
addition  to  tether  tension  control)  to  improve  system 
performance.  Additional  control  will  be  required  to  provide 
for  out-of-plane  damping  of  rigid  motions  and  vibration 
suppression . 

Because  most  of  the  useful  missions  are  carried  out  during 
the  station  keeping  phase  for  Shuttle/Platform  Tethered 
Subsatellite  systems,  a  review  of  the  various  tether  system 
control  laws  has  focused  mainly  on  the  deployment  and 
station  keeping  stages.  Retrieval  is  less  important  than 
the  first  two  for  tether  reflector  applications  where  it  may 
not  be  required  to  retrieve  the  tether  (except  possibly 
before  rapid  maneuvering) .  Although  some  nonlinear  control 
laws  were  proposed  (especially  those  which  include  thruster 
augmentation  of  the  tether  tension  control),  it  may  still  be 
difficult  to  implement  an  efficient  and  successful  retrieval 
for  the  tethered  reflector  system  and  further  study  is 
suggested . 

For  deployment,  the  active  tension  modulation  schemes 
proposed  by  Rupp  and  improved  by  subsequent  investigators 
are  mere  efficient  than  the  purely  passive  scheme  advocated 
by  Kane  and  Levinson.  A  tension  control  law.  based  on  a 
combination  of  exponential  and  uniform  rate  of  tether 
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length,  as  a  commanded  length  rate  can  be  used  both  for  in- 
plane  and  out-of-plane  motion  control  during  deployment. 

Out  of  the  various  active  control  laws  for  proposed  station 
keeping  of  orbiting  tethered  reflector  systems,  Rupp's 
control  law  is  the  most  basic  and  is  effective  in 
controlling  in-plane  motion,  but  not  adequate  for  out-of¬ 
plane  motion  control.  Control  laws  based  on  optimal  control 
theory  offer  the  greatest  potential  for  applications 
involving  proposed  orbiting  tethered  reflector  systems. 
Alternate  tension  modulation  optimal  control  laws  based  on 
both  in-plane  tether  swing  angle  and  vibrational  state 
information  can  result  in  further  improvements  as  compared 
with  Rupp's  control  law. 

For  out-of-plane  motion  control  during  station  keeping 
a  combination  of  tension  modulation  in  the  tether  plus  other 
forms  of  control  (such  as  the  use  of  thrusters)  will  be 
required . 

Finally  ,  a  preliminary  model  of  the  nonlinear  dynamics  of 
the  orbiting  tethered  antenna/reflector  system  has  been 
developed  based  on  Lagrangian  formulation.  The  resulting 
equations  are  highly  coupled  and  for  deployment  represent  a 
set  of  non-autonomous  differential  equations.  For  this 
model  the  shell  was  considered  to  be  rigid,  but  the  mass  and 
flexibility  of  the  tether  has  been  taken  in no  account. 

These  equations  will  be  used  in  the  next  phase  of  this 
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effort  to  simulate  deployment  dynamics  and  compare  the 
performance  using  different  control  strategies. 
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